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Abstract. We present optical studies of the physical and wind properties, plus CNO chemical abundances, of 25 09.5-B3 
Galactic supergiants. We employ non-LTE, line blanketed, extended model atmospheres, which provide a modest downward 
revision in the effective temperature scale of early B supergiants of up to l-2kK relative to previous non-blanketed results. 
The so-called 'bistability jump' at Bl (r cff ~ 21kK) from Lamers et al. is rather a more gradual trend (with large scatter) 
from Voo/Vesc ~3.4 for B0-0.5 supergiants above 24kK to v«,/v csc ~2.5 for B0. 7-1 supergiants with 20kK< T c g <24kK, and 
Voo/Vesc ~l-9 for Bl.5-3 supergiants below 20kK. This, in part, explains the break in observed UV spectral characteristics 
between B0.5 and B0.7 subtypes as discussed by Walborn et al. We compare derived (homogeneous) wind densities with recent 
results for Magellanic Cloud B supergiants and generally confirm theoretical expectations for stronger winds amongst Galactic 
supergiants. However, winds are substantially weaker than predictions from current radiatively driven wind theory, especially 
at mid-B subtypes, a problem which is exacerbated if winds are already clumped in the Her line forming region. In general, 
CNO elemental abundances reveal strongly processed material at the surface of Galactic B supergiants, with mean N/C and 
N/O abundances 10 and 5 times higher than the Solar value, respectively, with HD 2905 (BC0.7 la) indicating the lowest degree 
of processing in our sample, and HD 152236 (B1.5 Ia + ) the highest. 
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1. Introduction 

Early-type stars possess powerful winds that provide signifi- 
cant feedback to the local ISM (e.g. Smith 2005). Significant 
progress has recently been made towards an improved un- 
derstanding of the physical parameters and winds of mas- 
sive O stars and their immediate descendants, OBA super- 
giants. In particular, the availability of modern spectral synthe- 
sis codes allowing for line blanketing and stellar winds together 
with high-resolution instruments on 4-8m telescopes permits 
detailed studies of OB stars in the Magellanic Clouds (e.g. 
Crowther et al. 2002; Hillier et al. 2003; Trundle et al. 2004). A 
primary motive behind such studies is the investigation of how 
winds in the low-metallicity LMC and SMC differ from those 
in the Milky Way. 

For O stars, Herrero et al. (2003) and Repolust et al. (2004) 
provide suitable Galactic comparison stars, whilst one cur- 
rently has to rely on results from Kudritzki et al. (1999) for B 
supergiants, for which the temperature scale of McErlean et al. 
(1999) was adopted. This is one of the principal reasons behind 
the present study, namely a determination of the wind prop- 
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erties and temperature scale of B0-3 supergiants using con- 
temporary line blanketed techniques and optical spectroscopy. 
Although the range in spectral type is narrow, the effective 
temperatures of B0-3 supergiants span almost a factor of two, 
from 17 to 32kK (McErlean et al. 1999). The present optical 
study of ionization in early B supergiants complements the re- 
cent UV wind ionization investigation by Prinja et al. (2005). 
Further, there is now considerable qualitative evidence support- 
ing CNO-processed material in the atmospheres of massive 
stars (Herrero & Lennon 2004), although quantitative results 
for B supergiants - for comparison with recent evolutionary 
models allowing for rotation (Meynet & Maeder 2000) - are 
generally lacking. 

A second motive for the present study is an observational 
investigation of the so-called 'bistability jump' in the terminal 
velocities of early B supergiants (Lamers et al. 1995). Walborn 
& Nichols-Bohlin (1987) first remarked upon the UV morpho- 
logical differences between B0. 5 la and B0. 7 la supergiants, in 
the sense that Cn /U335 and Aim /11855-63 P Cygni profiles 
first appear at B0.7, whilst the Siiv /U393-1402 doublet be- 
comes narrower and stronger. In solely two cases, HD 115842 
(B0.5Ia) and HD 91943 (B0.7Ia) are these distinctions re- 
versed (Walborn et al. 1995). Pauldrach & Puis (1990) first 
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Table 1. Program Stars - observations are from either the lm JKT (Lennon et al. 1992), 2.5m INT or 4.2m WHT (northern) 
or 1.5m CTIO (southern). Two estimates of interstellar reddenings are given, based on intrinsic colours adopted from Schmidt- 
Kaler (1982) and from theoretical continuum fits to UV spectrophotometry and optical photometry (no spectrophotometry was 
available for HD 194279 or HD 148688). The majority of our targets are probable cluster/association members, whilst absolute 
magnitudes are adopted for the remainder, based on average values of Magellanic Cloud members from Fitzpatrick ( 1 99 1 ), Parker 
(1993) and Massey et al. (1995), except that HD 190603 is set at My - -7.5 mag. References to cluster/association distances are: 
(a) Brown et al. (1994); (b) Blaha & Humphreys (1989); (c) Baume et al. 1999; (d) Baume et al (2003); (e) Garmany & Stencel 
(1992); (f) Humphreys (1978). 



HD 


Name 


Sp Type V 


B-V 


E(B-V) 


Assoc. 


DM 


Ref 


My 


AM V 


Obs 






mag 


mag 


(B-V) 


UV 


mag 




mag 


mag 






30614 


a Cam 


09.5 la 4.29 


+0.02 


0.29 


0.32 




10.0 




-6.6 


±0.5 


WHT 


37128 


eOri 


B0 la 1.70 


-0.19 


0.05 


0.06 


Ori OB lb 


7.8 


a 


-6.3 


±0.5 


WHT/INT 


91969 




B0 la 6.52 


+0.00 


0.24 


0.26 


NGC3293 


12.2 


d 


-6.4 


±0.3 


CTIO 


94909 




B0 la 7.34 


+0.48 


0.72 


0.754: 




11.5 




-6.4 


±0.5 


CTIO 


122879 




B0 la 6.42 


+0.12 


0.36 


0.37 




11.7 




-6.4 


±0.5 


CTIO 


38771 


k Ori 


B0.5 la 2.05 


-0.18 


0.04 


0.07 


OriOBlc 


8.0 


a 


-6.1 


±0.5 


WHT/INT 


115842 




B0.5 la 6.02 


+0.29 


0.51 


0.57$ 




11.3 




-6.9 


±0.5 


CTIO 


152234 




B0.5Ia(Nwk) 5.44 


+0.20 


0.42 


0.45 


Sco OBI 


11.5 


c 


-7.4 


±0.5 


CTIO 


2905 


k Cas 


BC0.7Ia 4.16 


+0.14 


0.35 


0.39 


Cas OB 14 


10.2 


b 


-7.1 


±0.3 


WHT/INT 


91943 




B0.7Ia 6.73 


+0.07 


0.28 


0.22$ 


NGC 3293 


12.2 


d 


-6.3 


±0.3 


CTIO 


152235 




B0.7Ia(Nwk) 6.34 


+0.51 


0.72 


0.75$ 


Sco OBI 


11.5 


c 


-7.4 


±0.5 


CTIO 


154090 




B0.7Ia 4.87 


+0.26 


0.47 


0.49$ 




10.2 




-6.8 


±0.5 


CTIO 


13854 




Bllab 6.48 


+0.28 


0.47 


0.55 


Per OBI 


11.8 


c 


-6.8 


±0.5 


WHT/INT 


91316 


p Leo 


BlIab(Nstr) 3.85 


-0.14 


0.05 


0.06 




10.5 




-6.8 


±0.5 


CTIO 


148688 




Blla 5.33 


+0.33 


0.52 


0.52$ 




10.5 




-6.8 


±0.5 


CTIO 


14956 




B1.5Ia 7.19 


+0.72 


0.90 


0.95$ 


Per OBI 


11.8 


c 


-7.4 


±0.5 


WHT/INT 


152236 


C Sco 


B1.5Ia + 4.73 


+0.48 


0.66 


0.69 


Sco OBI 


11.5 


c 


-8.8 


±0.5 


CTIO 


190603 




B1.5Ia + 5.65 


+0.54 


0.72 


0.80$ 




10.9 




-7.5 


±1.0 


WHT/INT 


14143 




B2Ia 6.66 


+0.50 


0.67 


0.70 


ft Per 


11.8 


c 


-7.2 


±0.3 


WHT/INT 


14818 


10 Per 


B2 la 6.25 


+0.30 


0.47 


0.52 


Per OBI 


11.8 


c 


-7.0 


±0.5 


WHT/INT 


41117 


X 2 Ori 


B2Ia 4.63 


+0.28 


0.45 


0.48 


Gem OBI 


10.9 


b 


-7.6 


±0.3 


WHT 


194279 




B2Ia 7.05 


+ 1.02 


1.19 


1.26$ 


Cyg OB9 


10.4 


f 


-7.0 


±0.3 


WHT/INT 


198478 


55 Cyg 


B2.5Ia 4.86 


+0.42 


0.57 


0.59 


Cyg OB7 


9.5 


f 


-6.4 


±0.3 


WHT 


14134 




B3Ia 6.55 


+0.45 


0.58 


0.65$ 


h Per 


11.8 


c 


-7.1 


±0.3 


IKT/INT 


53138 


o 2 CMa 


B3Ia 3.01 


-0.08 


0.05 


0.07 




10.2 




-7.3 


±0.5 


IKT/INT 



$: UV spectrophotometry absent or poor quality, leading to uncertain reddening from this technique. 



described a bi-stability jump with reference to the extreme B 
supergiant P Cygni (B1.5 Ia + ), in the sense that a small change 
in physical properties for P Cyg can result in dramatic differ- 
ences in the outer wind, either ionized or recombined (Najarro 
et al. 1997; Crowther 1997). Amongst normal B supergiants, 
with much weaker winds, Lamers et al. (1995) demonstrated 
that the ratio Vco/v e sc drops steeply from about 2.6 for tempera- 
tures above 21kK (Bl subtype), to a value of 1.3 below 1 . 

The empirical B supergiant bistability jump of Lamers et al. 
has been widely used in synthetic spectral model calculations, 
of application to distant star forming galaxies (e.g. Rix et al. 
2004). Vink et al. (1999, 2000) predicted a jump in mass-loss 
around 25kK, due to the line acceleration of Fe m, such that 
the mass-loss increases by a factor of 5 between 27.5kK and 
22.5kK, although the change in velocity was not thoroughly 
explained. Here, we present studies of B supergiants to inves- 
tigate the Lamers et al. bistability jump, and test whether the 
Vink et al. predictions are consistent with observations. 



1 Prinja & Massa (1998) suggested a larger observational spread in 
Voo/Vesc amongst early B supergiants 



2. Observations 

We have selected 25 Galactic early-type supergiants in the 
spectral range 09.5-B3, 75% of which have established clus- 
ter/association membership, listed in Tabled In all but two 
cases, robust UV wind velocities are available from the liter- 
ature - otherwise average subtype values from Howarth et al. 
are adopted (HD 19456, HD 192479). Cluster or association 
distance estimates are available for the majority of the sam- 
ple, whilst average subtype absolute magnitudes are adopted 
for the remainder based on Magellanic Cloud stars (Fitzpatrick 
1991, Parker 1993, Massey et al. 1995). Many of the associ- 
ations are rather extended (e.g. Per OBI: Kendall et al. 1995) 
so in general we adopt uncertainties of +0.5 mag in distance 
modulus/absolute magnitude for field and association mem- 
bers. Reduced uncertainties of ±0.3 mag are restricted to clus- 
ter members (NGC 3293, h Per) and the more spatially con- 
centrated OB associations (Cyg OB7, Cyg OB9, Gem OBI, 
Cas OB 14). For the two Ori OBI members, we adopt distances 
from individual Ori OBI subgroups from Brown et al. (1994), 
lb for HD 37128 and Ic for HD 38771. The latter distance of 
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Table 2. Atomic models for early and mid B supergiants. For each ion F denotes full levels, S super levels, and T the number of 
bound-bound transitions considered. Most ions are common to all stars, except that high ionization stages are only included in 
early B supergiant models (shown in italics). 



Element I II III IV V 

FST FST FST FST FST 

H 30 20 435 

He 59 41 590 30 20 435 

C 53 30 323 54 29 268 18 13 76 

N 22 10 59 41 21 144 70 34 430 

O 75 18 450 63 22 444 45 25 182 

Mg 45 18 362 

Al 44 26 171 65 21 1452 

Si 62 23 365 45 25 172 12 8 26 

S 87 27 786 41 21 177 92 3 7 708 

Ca 12 7 28 

Fe 510 100 7501 607 65 5482 272 48 3113 182 46 1 781 



400pc is in preference to 220pc as derived from Hipparcos, 
given the difficulty in determining individual distance estimates 
along this line of sight (de Zeeuw et al. 1999). For HD 190603, 
we adopt an absolute magnitude of -7.5+1 .0 mag, intermediate 
between normal la supergiants and the hypergiant HD 152236. 
Finally, we consider the Hipparcos membership of HD 53138 
to Col 121 to be dubious since it has insignificant proper mo- 
tion compared to the cluster (de Zeeuw et al. 1999), and so 
instead adopt an average B3 la absolute magnitude. 

Interstellar reddenings (assuming R = Ay/E(B - V) = 
3.1 throughout) are determined using intrinsic colours from 
Schmidt-Kaler (1982) and independently from theoretical flux 
distributions determined in the present analysis, together with 
archival large aperture IUE spectrophotometry (except for 
HD 194279 and HD 148688). Although the former was 
adopted for distance modulus/absolute magnitude estimates, 
the latter provided a check on the reliability of intrinsic colours 
for early B supergiants. In general the combined UV/optical 
flux distribution suggested intrinsic (B-V)o values ~0.02 mag 
more negative than Schmidt-Kaler (1982) for our la sample. 

Optical observations originate from either the Isaac Newton 
Group, La Palma (northern sample) or CTIO, Chile (south- 
ern sample). The northern targets were originally observed 
in the blue visual with the JKT (Lennon et al. 1992) with 
the Richardson-Brealey Spectrograph with the R1200B grat- 
ing at a resolution of 0.8A and S/N~150. The majority were 
re-observed with the ISIS spectrograph at the WHT at a simi- 
lar resolution (Smartt, priv. comm.) though higher S/N(>300) 
using a setup identical to that used in Smartt et al. (2001). Ha 
observations were obtained at the INT using the IDS spectro- 
graph by DJL, with a spectral resolution of 1 .2A, and S/N~200. 
In three cases, HD 30614, HD 41117 and HD 198478, WHT 
ISIS observations were obtained by PAC during 2002 October 
using a dichroic and 1.2arcsec wide slit, providing simultane- 
ous blue and Ha spectroscopy with the 1200B/R gratings in 
first order, centred at i4400 and ,16650, with S/N>300. 

The southern sample was observed at the CTIO 1.5m dur- 
ing May-June 2003 via the SMARTS service program using 



the Cassegrain spectrograph with Loral IK (1200x 800 pixel) 
CCD. A standard reduction, bias subtraction, flat-fielding, and 
extraction were carried out using IRAF. Blue and red spec- 
troscopy of each target was obtained using the 47 grating, pro- 
viding a spectral range, and (3 pixel) spectral resolution of 
660A, 1.6 A (second order, blue) and 1320 A 3. OA (first order, 
red) respectively, as measured from HeAr and Ne arc lines. The 
extracted blue and red continuum S/N ranges from 150 to 200, 
and 200 to 300 respectively. 

In order to assess the robustness of our analyses based 
on intermediate resolution observations covering a restricted 
wavelength range, HD 1 15842 was observed for us by Dr C.J. 
Evans (ESO programme 073.D-0234) at the MPG/ESO 2.2m 
on 7 July 2004 with the Fibre-fed, Extended Range Echelle 
Spectrograph (FEROS). A two arcsec aperture provides a re- 
solving power of R =48, 000 covering a spectral range of 3500- 
9200 A using an EEV 2kx4k CCD. A single exposure of 240s 
provided a continuum S/N of 300 at /14600 in the pipeline re- 
duced dataset (A. Kaufer, priv. comm.). 

Finally, (at least) two of the present sample definitely pos- 
sess significant Ha variability, HD 37128 (B0 la, Prinja et al. 
2004) and HD 14134 (B3 la, Morel et al. 2004). For HD 37128 
we investigate the effect of variability on the derived mass-loss 
rates using Ha observations provided by Dr R.K. Prinja. 



3. Analysis 

3.1. Methods 

Recent developments in model-atmosphere studies of hot stars 
routinely permit sophisticated metal-line blanketing in non- 
LTE to be considered, under either plane-parallel geometry 
(e.g., TLUSTY: Hubeny & Lanz 1995; Lanz & Hubeny 2003) 
or spherical geometry (e.g., CMFGEN: Hillier & Miller 1998). 
The former technique provides a more comprehensive treat- 
ment of photospheric metal-line blanketing, while the latter of- 
fers insights into wind properties. 



4 



Paul A. Crowther et al.: Early B supergiants 




Temperature (kK) 



Fig. 1. logg - T e ff calibration from Trundle et al. (2004) and 
Trundle & Lennon (2005) for SMC B supergiants (open trian- 
gles) plus HD 30614 from Repolust et al. (2005, open circle). 



In the present application, we have utilized CMFGEN, the 
current version of which is discussed by Hillier et al. (2003). 
The code solves the radiative transfer equation in the co- 
moving frame, under the additional constraints of radiative and 
statistical equilibrium. The model used here is similar to that 
applied by Evans et al. (2004b). In total, 26 ions of H, He, C, 
N, O, Mg, Al, Si, S, Ca, Fe are included for early subtypes, 
comprising a total of 2,636 individual levels (grouped into 766 
superlevels), with a full array of 25,960 bound-bound transi- 
tions. Higher ionization stages are omitted for B2.5-3 super- 
giants as indicated in Table |2] 

The temperature structure is determined by radiative equi- 
librium. CMFGEN does not solve the momentum equation, so 
that a density or velocity structure is required. For the super- 
sonic part, the velocity is parameterized with a classical y6-type 
law, with the exponent in the range 0=1-3 selected on the ba- 
sis of the Ha profile. This is connected to a hydrostatic density 
structure at depth, such that the velocity and velocity gradient 
match at the interface. The subsonic velocity structure is set 
by the corresponding H-He TLUSTY model which were cal- 
culated by the authors for each CMFGEN model atmosphere 
calculation. With regard to CMFGEN models in which the 
TLUSTY structure at depth is neglected, stellar temperatures 
are ~lkK lower and mass-loss rates 0.1-0.2 dex higher. 

We initially estimated surface gravities at individual tem- 
peratures from a 2nd order fit to results from Trundle et al. 
(2004) and Trundle & Lennon (2005) for la B supergiants, 
plus HD 30614 (09.51a) from Repolust et al. (2004). This is 
presented in Fig. ^ an d should lead to gravities accurate to 
+0.2 dex. Subsequently, we have adjusted surface gravities in 
TLUSTY models until Hy is reproduced. We anticipate accu- 
racies of 0.1-0. 15 dex for all targets except perhaps the hyper- 
giants HD 152236 and HD 190603 which have significant wind 
contamination at Hy. 

We have assumed a depth-independent Doppler profile for 
all lines when solving for the atmospheric structure in the co- 
moving frame, while in the final calculation of the emergent 



spectrum in the observer's frame, we have adopted a radially 
dependent turbulence 2 , which reflects the effect of shocks due 
to wind instabilities. For the present sample we initially assume 
£ = 20 kms -1 , with other values in the range 10^-0 kms -1 con- 
sidered if consistent fits are not obtained for the helium and 
silicon lines. Incoherent electron scattering and Stark broaden- 
ing for hydrogen and helium lines are adopted. We convolve 
our synthetic spectrum with a rotational broadening profile. 
Because of the intermediate dispersion of most observations, 
instrumental effects prevent a reliable determination of v sin i 
from the present dataset, such that we adopt the UV derived 
rotational rates from Howarth et al. (1997). 

Terminal velocities, Vco, were uniformly determined from 
Vbiack measurements from the literature (Howarth et al. 1997), 
and are listed in Table |4] Alternatively, velocities may be 
derived via the Sobolev with Exact Integral (SEI) technique 
(Lamers et al. 1987; Haser 1995). Although individual mea- 
surements differ by as much as 20%, v 00 (Haser)/v 00 (Howarth 
et al.)=1.02 for the 19 stars in common to both samples. SEI 
modelling is often preferred to Vbiack measurements of wind ve- 
locity, although has its own modelling uncertainties. For ex- 
ample, Haser (1995) and Evans et al. (2004a) both used SEI 
modelling of UV lines in HD 152236, deriving v,x,=500 km s _1 
and 450 kms -1 , respectively, versus 390 kms _I adopted here 
from Howarth et al. (1997). 

We varied the mass-loss rate and the velocity law (as char- 
acterized by the exponent /3) until the Ha profile was best re- 
produced in terms of intensity and morphology. Next, we ad- 
justed the stellar temperature to match the intensities of the 
He ii /14542 and Hei A4471 lines for O9.5-B0 stars, Si IV 
,14089 and Si III A4553-68-75 lines for B0.5-B2 stars, and Si II 
,14128-31 and Si III /14553-68-75 lines for B2.5-3 subtypes. 
The accuracy with which we can derive stellar temperatures is 
+ 1000K. Radii and luminosities follow from our derived ab- 
solute magnitude, whilst mass-loss rates are obtained from a 
spectral fit to the Ho- profile. With absolute magnitudes fixed 
at our adopted values, formal uncertainties of 10% in dM/dt 
follow for stars with strong winds, increasing to 30% for weak 
lined stars. In general, the greatest uncertainty relates to the /3 
law, although the Ha profile does provide some constraints on 
the p exponent. 

He/H abundances are relatively difficult to constrain in B 
supergiants since different choices of Hei lines yield differ- 
ent He abundances, in part due to the well known singlet- 
triplet problem (Lennon et al. 1991; McErlean et al. 1998). 
Consequently, we adopt a uniform He/H = 0.2 by number for 
all our targets. This was selected on the basis of the average 
N/C enrichment derived relative to Solar values (Sect. 14. 2> . 
with respect to recent interior evolutionary models for rotating 
massive stars (Meynet & Maeder 2000). Fortunately, an error 
in He/H will only weakly affect the derived stellar temperature 
and/or mass-loss rates. Indeed, Repolust et al. (2005) have re- 
cently obtained a typical He/H abundance of 0.2 amongst their 
sample of Galactic OB supergiants. 



2 In CMFGEN this adopted value is technically £ rain , the microtur- 
bulence at the base of the wind, which increases with radius to £ max = 
50 kms~' at v^. 



Paul A. Crowther et al.: Early B supergiants 



5 



Table 3. Derived non-LTE stellar and wind parameters. Two values are listed for the mass-loss rate of HD 37128 since we fit the 
minimum and maximum Ho- profiles of Prinja et al. (2004). 



HD 


Name 


Sp Type 


Teg 


logg 


loggeff 


log L/L G 


R, 


Voo 


Voc,/Ve sc 


Mx 10 6 


P 


v sin i 




My 


MvocR 05 








kK 


cgs 


cgs 




Re 


kms~' 




M /yrkms 1 


kms~' 


kms~' 


mag 


cgs 


30614 


a Cam 


09.5 la 


29.0 


3.0 


2.75 


5.63 


26.0 


1560 


3.40 


5.0 


1.5 


129 


30 


-6.6 


29.40 


37128 


e Ori 


B0 la 


27.0 


2.9 


2.85 


5.44 


24.0 


1910 


4.76 


(max) 2.5 
(min) 2.0 


1.5 


91 


12.5 


-6.3 


29.17 
29.07 


91969 




B0 la 


27.5 


2.95 


2.75 


5.52 


25.3 


1470 


3.33 


1.0 


1.5 


79 


7.5 


-6.4 


28.67 


94909 




B0 la 


21.0 


2.9 


2.7 


5.49 


25.5 


1050 


2.55 


2.0 


1.5 


100£ 


10 


-6.4 


28.62 


122879 




B0 la 


28.0 


2.95 


2.7 


5.52 


24.4 


1620 


3.81 


3.0 


1.5 


92 


15 


-6.4 


29.18 


38771 


k Ori 


B0.5 la 


26.5 


2.9 


2.7 


5.35 


22.2 


1525 


3.84 


0.9 


1.5 


83 


12.5 


-6.1 


28.61 


115842 




B0.5 la 


25.5 


2.85 


2.65 


5.65 


34.2 


1180 


3.86 


2.0 


1.5 


84 


10 


-6.9 


28.94 


152234 




BO 5 la (N wkl 


26.0 


2.85 


2.6 


5.87 


42.4 


1450 


2.85 


2.7 


1.5 


100$ 


10 


-7.4 


29.21 


2905 


k Cas 


BC0 7 la 


21.5 


2.6 


2.45 


5.52 


41.4 


1105 


2.79 


2.0 


2.0 


91 


20 


-7.1 


28.95 


91943 




B0.7Ia 


24.5 


2.8 


2.6 


5.35 


26.3 


1470 


3.72 


0.75 


1.2 


79 


10 


-6.3 


28.55 


152235 




BO 7 la IN wkl 


23.0 


2.65 


2.45 


5.76 


47.1 


850 


1.98 


1.25 


1.5 


81 


10 
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p. Our study requires higher v sin i values for HD 94909 and HD 152234 than 68 and 76 km s 1 obtained by Howarth et al. (1997), respectively, 
whilst HD 194279 and HD 14956 were not included in their sample. 



CNO abundances were varied to fit the relevant optical 
lines, with a typical accuracy of ~0.2 dex. Solar abundances 
were taken from Grevesse & Sauval (1998) except oxygen, 
which was set at log(0/H)+12 = 8.66 (Asplund et al. 2004). 
Abundances of other metallic elements were fixed at Solar 
values. The primary diagnostic lines used for determination 
of the nitrogen abundance were N II ,13995 and ,14601-43, 
plus N III ,14097 blend for the earliest subtypes where N II is 
weak/absent. Oxygen and carbon abundances were also deter- 
mined from consideration of the optical lines; O II ,14069-4092 
and ,14590-96 provided useful constraints on the oxygen abun- 
dance, whilst C II ,14267 together with C III /14647-51 at early 
subtypes (blended with O II) were used for carbon. We antic- 
ipate accuracies in derived CNO abundances to be no greater 
than ~0.15 dex (O) to ~0.3 dex (C and N), where differences 
are due to the number of diagnostics available. 

4. Spectroscopic Results 

Derived stellar, wind and chemical properties of our sample 
of OB supergiants based on the intermediate dispersion spec- 
troscopy are provided in Table[5] We include surface and effec- 
tive gravities, corrected for the Eddington parameter, Y e , which 
is the ratio of the radiative acceleration from electron scattering 
to the stellar gravity. 



4.1. Physical and wind properties 

Figure |2ja-c) compares blue optical spectroscopy for selected 
09.5-B3Ia supergiants with synthetic spectral fits. In general, 
the quality of spectroscopic fits is excellent, not only for our se- 
lected He or Si diagnostic lines, but also for classification lines 
such as He i ,1447 1 vs Mg n ,1448 1 , plus weak features of other 
non-CNO elements, e.g. Aim /14529 and Siv ,14485-4504. 
The major differences reflect weak metal transitions, which are 
missing from our selected model atoms, with no effect upon 
the derived physical properties. Apparent discrepancies in the 
region of ,14430 are in fact due to the broad Diffuse Interstellar 
Band (DIB) at that wavelength, whose strength scales with in- 
terstellar reddening. The only clear discrepancies occurs for 
HD 30614 (09.51a) and HD 37128 (B0 la) at Hen ,14686 and 
Si iv ,14088-4116. 

In Fig. [3] we present spectral fits to all our program stars in 
the vicinity of Ha, with model profiles degraded to the spectral 
resolution of the individual observations. Overall, agreement is 
excellent, except that the blue P Cygni absorption component is 
imperfectly reproduced for the Bl-2 supergiants. Note that the 
Cn ,1,16578-82 doublet is well matched in mid B supergiants 
in most cases. As already indicated, HD 37128 (B0 la) is ob- 
served to exhibit large scale Ha variability (Prinja et al. 2004) 
so we have attempted to reproduce the minimum and maximum 
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Fig. 2. (a) Synthetic fit (red in electronic version, dotted in paper version) to /L13975^4250 spectral range for representative early 
B supergiants (solid black). 



emission cases to quantify the range in derived mass-loss rate, Table 0] compares our derived temperatures for B super- 
namely 2.0-2.5 xl0~ 6 M yr . giants with those derived previously by McErlean et al. (1999), 
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Fig. 2. (b) Synthetic fit to AA42 15^4485 spectral range for representative early B supergiants. 



using hydrostatic, plane-parallel model atmospheres. As with line-blanketed, spherically extended model atmospheres indi- 
O stars (e.g. Crowther et al. 2002, Repoulst et al. 2004), current cate a downward revision in the temperature scale, by typically 
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l-2kK Indeed, stellar temperatures obtained with CMFGEN indicated by Urbaneja (2004) who has used FASTWIND to 
seem to be agree well with FASTWIND (Puis et al. 2005), as study HD 37128, HD 38771, and HD 14956 revealing tem- 




Fig. 3. Spectral fits (red in electronic version, dotted in paper version) to Ha and C n -16578—82 for 09.5-B3 supergiants. 



peratures consistent with our results to ~0.5kK. Given this 
close agreement, we compare bolometric corrections (BCs) as 
a function of temperature from the present sample, supple- 
mented by SMC B supergiant results from Trundle et al. (2004) 
and Trundle & Lennon (2005) in Fig@] A fit to the combined 
dataset indicates 

BC = 20.15 -5.13 log r eff 

with a standard deviation of ~0.05 mag, which closely matches 
the recent result for O stars by Martins et al. (2005) at ~30kK. 

Spectroscopically estimated masses span a wide range, 
from 8^M Q (HD 198478) to 76+£*M (HD 152236), based on 
gravities accurate to logg+0.15 dex. In general, spectroscopic 
masses are lower than estimates of the evolutionary masses, 
indicating another case of a 'mass discrepancy', as already es- 
tablished by Trundle et al. (2004) for SMC B supergiants. For 
example, we estimate an initial (current) evolutionary mass of 
2OM (18 M ) for HD 198478. In order to resolve the spectro- 
scopic and evolutionary masses, we would require \ogg ~2.5, 
versus 2. 15+0. 15 obtained for this star. 

In Table 13 we compare the present results with those from 
previous studies, namely (i) Barlow & Cohen (1977) based 
on the mid-IR free-free excess technique of Wright & Barlow 
(1975) and adopted temperature scale for B supergiants, and 
(ii) Kudritzki et al. (1999) who adopted the temperature scale 
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Fig. 4. Bolometric correction (BC) as a function of temperature 
for Galactic (filled circles, this study) and SMC (open triangles, 
Trundle et al. 2004; Trundle & Lennon 2005) B supergiants. 
The fit to the combined sample (solid line) connects well to the 
recent Martins et al. (2005) relation for O stars (dotted line). 



10 Paul A. Crowther et al.: Early B supergiants 



of McErlean et al. (1999), based on hydrostatic, plane-parallel 
model atmospheres and derived Ha mass-loss rates with the un- 
blanketed version of FASTWIND (Santolaya-Rey et al. 1997); 
(iii) Urbaneja (2004) based on FASTWIND (Puis et al. 2005) 
model atmospheres. 

Table 4. Comparison between our derived non-LTE line blan- 
keted temperature scale (in kK) of Galactic early B la and lab 
supergiants and the unblanketed temperature scale of McErlean 
et al. (1999), plus recent line blanketed Magellanic Cloud re- 
sults from Trundle et al. (2004), Trundle & Lennon (2005). For 
each study the range in temperatures and number of stars stud- 
ied is indicated. Values for hypergiants are flagged with *, and 
are excluded from our current "mean" temperature scale for la 
supergiants. 



Subtype McErlean et al. Trundle et al. This work T eff 
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Kudritzki et al. (1999) adopted the Hipparcos distance to 
HD 38771, whilst we adhere its the membership of Ori OBlc. 
Putting HD 38771 and the Ha variable HD 37128 aside, con- 
sistency between Kudritzki et al. and the present sample is 
excellent for the other B0-1 supergiants in common, with 
Mna/M^gi) ~ 1+0.1. In contrast, the present results are sys- 
tematically higher for B2-3 supergiants, with M Hq ./M K 99 ~ 3, 
except for HD 53138 for which substantially different abso- 
lute magnitudes were adopted in the two studies. This result, 
together with the good consistency between our results using 
CMFGEN and FASTWIND (Urbaneja 2004) suggests prob- 
lems with the application of the analytical mass-loss technique 
(originally devised for O stars by Puis et al. 1996) to B2-3 su- 
pergiants by Kudritzki et al. (1999). 

For those 10 stars in common with Barlow & Cohen (1977), 
if we adjust their results to our adopted wind velocities and 
distances, we find ratios M H a./M m id-iR in the range 0.7-2.3 
with an average of 1.4, excluding HD 53138 for which the 
mid-IR continuum exceeds the Ha result tenfold. Kudritzki et 
al. (1999) also comment on the peculiar mid-IR brightness of 
HD 53138 with regard to optical and radio mass-loss indices. 
Overall, agreement is satisfactory between the two techniques. 
This is of particular interest, given clumping, if present, is pre- 
dicted to affect these measurements in different ways - Ha 
samples the highest density region immediately above the pho- 
tosphere, whilst the mid-IR continuum samples free-free emis- 
sion , superimposed upon the stellar photosphere at larger radii 
(Runacres & Owocki 2002). Spectroscopic studies in which 



Table 5. Comparison between present results for selected B 
supergiants and those determined from (a) mid-IR free-free 
excess method of Barlow & Cohen (1977, BC77), and (b) 
Kudritzki et al. (1999, K99) following the Ha fitting technique 
of Santolaya-Rey et al. (1997) plus the plane-parallel hydrostat- 
ically derived temperature scale from McErlean et al. (1999). 
For HD 30614 we compare the present results with both Puis et 
al. (1996, P96) and the recent line blanketed study of Repolust 
et al. (2004, R04), whilst we also include line blanketed re- 
sults from Urbaneja (2004, U04) for HD 37 1 28 and HD 3877 1 . 
Wind velocities in parenthesis relate to adopted values. 



HD Sp Type T eff log(L/L G ) M x 10 6 M v Ref 

kK km/s M Q yr~' mag 



clumping is considered generally produce superior spectral fits 
to far-UV lines (e.g. Evans et al. 2004b). Consequently, a com- 
parison with Barlow & Cohen (1977) suggests comparable 
clumping factors in the Ha and mid-IR continuum forming re- 
gions. 



30614 09.5 la 30.0 5.79 1550 5.2 -7.0 P96 

29.0 5.83 1550 6.0 -7.0 R04 

29.0 5.63 1560 5.0 -6.6 This work 

37128 BOIa 25.0 5.62 2300 1.7 -6.8 BC77 

28.5 5.86 1600 2.4 -7.0 K99 

27.5 5.69 1600 2.5 -6.9 U04 

27.0 5.44 1910 2-2.5 -6.3 This work 
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26.5 4.94 1350 0.48 -4.9 U04 
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22.0 5.47 1110 0.35 -6.8 This work 
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19.0 5.66 (1000) 2.9 -7.5 BC77 

19.5 5.63 485 2.5 -7.5 This work 

18.0 5.46 (1000) 1.2 -7.1 BC77 

20.0 5.51 650. 0.30 -6.95 K99 

18.0 5.42 645 1.05 -7.2 This work 

20.0 5.47 650. 0.25 -6.9 K99 

18.5 5.35 565 0.55 -7.0 This work 

18.0 5.58 (1000) 1.2 -7.4 BC77 

19.5 5.70 500 0.85 -7.5 K99 

19.0 5.65 510 0.9 -7.6 This work 

15.0 5.16 (580) 0.56 -6.8 BC77 

16.5 5.03 470 0.23 -6.4 This work 

15.0 5.24 (580) 0.71 -7.0 BC77 

18.0 5.48 465 0.15 -7.1 K99 

16.0 5.28 465 0.52 -7.1 This work 

15.0 5.28 580 1.7 -7.1 BC77 

18.5 5.22 620 0.095 -6.4 K99 

15.5 5.34 865 0.36 -7.3 This work 



2905 

13854 Bllab 

91316 Bllab 

14956 B 1.5 la 

190603 B1.5Ia + 
14143 



B2Ia 



14818 
41117 



B2Ia 



B2Ia 



198478 B2.5Ia 
14134 B3Ia 

53138 B3Ia 



Paul A. Crowther et 



5.0 



- [N/C] # 


Z=0.03 V=300 km/s ■ 


r\ y [N/C]=+1.3 

O / 








' \ ^r1% • 




; 1/ 1 €/. 


- 


; s ~> T ° 


- - - J : 


V E 5m/ X [N/C] ^°- B 


O [N/C]<41- 

\m [n/c>+i 


\/SOM " ; ; 'l5M@ 


4.6 4.4 4.3 


4.0 j.s s.e 


log Teff 





a.g - 



■ [N/O] • Z = 0.0S V=300 km/s - 

e.o - 




log Teff 



Fig. 5. (Upper panel) Observed H-R diagram for program 
stars, coded by degree of N/C enrichment, together with tracks 
from Solar metallicity, rotating (300 km/s) evolutionary mod- 
els (dotted lines) from Meynet & Maeder (2000), for which 
predicted N/C enrichments are indicated at T e g =22,000K. The 
error bar for HD 91943 is indicated (±0.3 mag uncertainty in 
distance modulus); (Lower panel) As above for N/O. The er- 
ror bar for HD 190603 is indicated (+1.0 mag uncertainty in 
distance modulus). 

4.2. Elemental abundances 

Elemental abundances are presented in Table [6] These repre- 
sent the first attempt to quantify metal abundances in a sub- 
stantial sample of Galactic early and mid B supergiants. With 
respect to current Solar values (Asplund et al. 2005), car- 
bon is depleted by 0.0-0.8 dex, nitrogen is enriched by 0.1- 
1 dex, with oxygen typically mildly depleted (up to 0.5dex). 
Consequently, the mean [N/C] and [N/O] ratios for our sample, 
+ 1.07 and +0.74 dex respectively, indicate evidence for partial 
(rotationally mixed?) CNO processing at their surfaces. 

These conclusions support the relatively few quantitative 
abundance estimates for B supergiants from the literature. 
Previous studies typically relied upon differential abundance 
analyses (e.g. Smartt et al. 1997) using blanketed LTE model 
atmospheres (Kurucz 1991), supplemented by non-LTE plane 
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Table 6. Individual and mean CNO abundances for Galactic B 
supergiants (given as log N{X)/N(H) + 12), where X represents 
the appropriate element, versus Solar abundances (Asplund et 
al. 2005), to the nearest 0.05 dex. For all stars He/H=0.2 by 
number was adopted. In addition, we show the mean non-LTE 
abundances obtained from two Bl-2Ib supergiants from Gies 
& Lambert (1992, GL92), four B2Ia supergiant results from 
Lennon (1994, L94), plus luminosity class Il-Ia Galactic A su- 
pergiants from Venn (1995, V95). 
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parallel, hydrostatic model atmospheres (e.g. Lennon et al. 
1991, Gies & Lambert 1992). We include results for four mor- 
phologically normal B2Ia supergiants (including HD 14818 
and HD 14143) from Lennon (1994), plus two Bl-2Iab/b su- 
pergiants (including HD 91316) from Gies & Lambert (1992) 
in Table [6] For stars in common between the studies, Gies & 
Lambert (1992) obtained [N/C]=+l.ll dex (we derive +1.45) 
for HD 9 1 3 1 6 (B 1 lab) and [N/C] =+ 1 . 1 3 dex (versus +0.98 dex 
here) for HD 198478 (B3Ia), whilst Lennon (1994) obtained 
[N/C] =+1.6 dex (we derive +1.28 dex) and [N/O] =+0.58 (ver- 
sus +0.78 dex) for HD 14818 (B2Ia), and similar differences 
for HD 14143 (B2Ia). Venn (1995) obtained rather lower ni- 
trogen enrichment in a sample of typically lower mass A su- 
pergiants and bright giants (see Table|6j. 

In Fig.|5]we present the observed H-R diagram for our pro- 
gram stars for [N/C] and [N/O] coded by the degree of enrich- 
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ment, together with initially rotating (300 km/s) evolutionary 
models for Solar metallicity from Meynet & Maeder (2000). In 
general, greater nitrogen enrichment is observed in lower tem- 
perature stars, whilst greater enrichment is predicted at higher 
luminosity. Note that the lower luminosity mid-B supergiants 
appear to show greater processing of oxygen than carbon from 
Fig.0 

For each of the 20, 25 and 40M o tracks we present the 
predicted enrichment at a temperature of 22kK. On average, 
the predicted [N/C] and [N/O] ratios of +0.8-1.3 and +0.7-1.2 
match observations rather well, with the highest luminosity star 
HD 152236 revealing the greatest nitrogen enrichment, also as 
predicted. 

Within our sample, all are morphologically normal B su- 
pergiants, except that k Cas has a BC subtype (Walborn 1972), 
suggesting a lower degree of CNO processing, In fact, HD 2905 
does have the least enriched N abundance of our sample (two 
N wk supergiants are also only moderately enriched), but does 
appear to be moderately depleted in C, with no significant pro- 
cessing in O. Consequently, HD 2905 does indeed possess the 
lowest [N/O] and [N/C] of our sample, confirming its morpho- 
logical difference, potentially due to a low initial rotation ve- 
locity. 

Finally, Urbaneja (2004) has used FASTWIND to derive 
physical properties and metal abundances of several Galactic 
early B supergiants, including HD 37128, HD 38771 and 
HD 14956. For the B0-0.5 la supergiants CNO elemental abun- 
dances are consistent at the 0.2 dex level between the two ap- 
proaches, except that Urbaneja (2004) obtains a 0.25 dex lower 
N abundance for HD 37128 and a 0.3 dex lower C abundance 
for HD 3877 1 . For HD 14956 (Bl .5 la), the agreement is poorer 
- Urbaneja (2004) obtains a 0.55 dex lower C abundance, 0.15 
dex lower N abundance and 0.4 dex higher O abundance. Such 
differences should be taken into account regarding the absolute 
uncertainty in abundances resulting from such studies. 



5. Wind properties 

Recently, Trundle et al. (2004) and Trundle & Lennon 
(2005) have analysed a sample of SMC B supergiants using 
FASTWIND, based upon high and intermediate dispersion op- 
tical spectroscopy. Before we compare the wind properties 
of Galactic and SMC B supergiants we wish to investigate 
whether our use of intermediate dispersion datasets introduces 
systematic effects on derived stellar properties. 

5.1. Analysis of high resolution observations of 
HD 115842 

HD 1 15842 was observed with both the CTIO 1 .5m/Cassegrain 
and ESO 2.2m/FEROS, allowing us to compare the derived 
properties of this star from the two datasets. In fact, we 
were able to match our intermediate dispersion model for 
HD 115842 to the FEROS echelle dataset, except that the tur- 
bulence was increased from v f = 10 kms -1 to 25 kms -1 . We 
compare selected regions of the FEROS dataset and synthetic 
spectrum in Fig.|6] including regions not covered by the CTIO 
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Fig. 6. Spectral comparison between FEROS spectroscopy of 
HD 115842 (solid black) and synthetic spectroscopy (red in 
electronic version, dotted in paper version), with CTIO spec- 
troscopy additionally included, offset by -0.2 continuum units. 
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dataset. Overall, the agreement is excellent, such that we are 
confident in the use of intermediate dispersion for our analysis. 

5.2. Wind momenta for Galactic versus SMC B 
Supergiants 

We present the wind momenta of our Galactic sample versus 
luminosity in Fig. separated between early B (>23kK) and 
mid B (<23kK) subtypes, together with current predictions for 
early and mid subtypes from Vink et al. (2000). In addition, we 
include recent Magellanic Cloud results of Evans et al. (2004b), 
Trundle et al. (2004) and Trundle & Lennon (2005). The effect 
of an uncertainty in distance modulus of ±0.5 mag is indicated 
for HD 38771. This typical uncertainty illustrates why we do 
not attempt a calibration of the wind momentum luminosity 
relationship at this time. 

For B0-0.5 subtypes, Galactic wind momenta lie close to 
the Vink et al. (2000) predictions for the majority of cases. 
In three cases, the derived wind momenta fall 0.5 dex below 
predictions, namely HD 91943, 115842 and 152234. Due to 
this spread, there is no clear distinction between the observed 
wind momenta of early B supergiants in the SMC and Milky 
Way (compare triangles to circles), although this could be due 
to small number statistics. If winds of early B supergiants are 
clumped in the Ha line forming region, we would need to shift 
the observed wind momenta to smaller values. 

For B0. 7-3 subtypes, the measured mass-loss rates now lie 
typically ~0.5 dex below the calibration of Vink et al., with 
the exception of HD 2905 (BC0.7 la). As above, if Ha derived 
mass-loss rates of OB stars need to be corrected for clumped 
winds, as suggested by recent observational evidence (Massa 
et al. 2003; Evans et al. 2004b), the difference between theory 
and observation would be exacerbated since the observed data 
points would move to lower values. In contrast with the case 
for the earliest B supergiants, we now see a clear separation 
between Galactic (circles) and SMC (triangles) B0. 7-3 super- 
giants at approximately the predicted offset. 

5.3. Bistability jump 

Let us now return to the question of the bistability jump. First, 
having revised the temperature scale for early B supergiants, 
let us investigate the sharp ~21kK bistability jump of Lamers 
et al. (1995) based on our present results. 

In the upper panel of Fig.[8]we repeat the original Lamers 
et al. (1995) 09.5-B3 supergiant sample, which illustrates the 
sharp 21kKjump. This clear step is due to Lamers et al. adopt- 
ing a single temperature of 20.8kK for all their seven Bl su- 
pergiants. In the lower panel of Fig. [3] we present our current 
results, where we have omitted HD 14956 and HD 194279 for 
which no reliable wind velocities are known. This sample in- 
cludes several sources in the original Lamers et al. study, to- 
gether with values from Evans et al. (2004a), based on mea- 
sured HST/STIS wind velocities plus model atmosphere results 
from Evans et al. (2004b), Trundle et al. (2004), Trundle & 
Lennon (2005). 



HO-0.5 la (>S3kK) 



Gal # 
LMC ■ 

SMC * 
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lag L/Lq 



6.0 8.3 



Fig. 7. Comparison between the Galactic B supergiant wind 
momenta presented here, separated between 'early' (>23kK) 
and 'mid' (<23kK) subtypes, together with other recent studies 
based on CMFGEN (filled, Evans et al. 2004b) or FASTWIND 
(open, Trundle et al. 2004; Trundle & Lennon 2005; Repolust 
et al. 2004; Herrero et al. 2002), plus radiatively driven wind 
theory predictions by Vink et al. (2000, 2001) for early and 
mid B supergiants. The effect of an uncertainty of +0.5 mag in 
distance modulus is indicated for HD 38771 (B0. 5 la). 



Error bars allow for uncertainties in distance, surface grav- 
ity and temperature. Since distances to Galactic early-type su- 
pergiants are notoriously imprecise (witness the Hipparcos ver- 
sus Ori OBI distance estimates to HD 38771), we provide error 
bars allowing for uncertainties in distance and surface gravity. 
Of these, the latter dominates, even for cases with particularly 
poorly constrained distances, such as HD 190603. 

In contrast with the earlier study, there is a gradual down- 
ward trend of v^/v^c with temperature, albeit with a large 
scatter. For T e ff > 24kK (approximately B0. 5 la and ear- 
lier), Voo/v esc ~ 3.4, for 20kK< T e ff <24kK (approximately 
B0. 7-1 la) Voo/v esc ~ 2.5, and for T e ff <20kK (approximately 
B 1.5 la and later) Voo/v eS c ~ 1-9- This reveals that the Bl 
'jump' is misleading in the context of normal B supergiants. 
Prinja & Massa (1998) came to similar conclusions based on 
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Fig. 8. The ratio Vco/vesc as a function of effective tempera- 
ture for (upper panel:) Galactic B supergiants from Lamers et 
al. (1995) together with the Kudritzki & Puis (2000) scalings 
(dotted lines); (lower panel): Galactic B supergiants from the 
present study (filled circles), together with results from Evans 
et al. (2004a, open triangles) for SMC B supergiants. Error bars 
relate to uncertainties in temperature (+lkK), gravity (+0.15 
to 0.2 dex) and absolute magnitude (+0.3 to 1 mag). We have 
omitted the two stars (HD 14956, HD 194279) for which reli- 
able wind velocities are not known. 



a larger B supergiant sample, albeit with an adopted subtype- 
temperature calibration. Amongst the present sample, the hy- 
pergiant HD 152236 (B1.5Ia + ) has the lowest ratio (1.0+°-f) 
whilst HD 37128 (B0 la) has the highest ratio of all (4.8+*°).' 

Vink et al. (1999) predicted a dramatic increase in wind 
densities of early B supergiants below ~25kK due to Fe, the 
dominant line driving ion, shifting to lower ionization stage. 
This temperature is close to the boundary between B0. 5 and 
B0. 7 subtypes where changes in v^jv^ do occur, although the 
original Lamers et al. (1995) 21kK boundary is an equally valid 
Voo/vesc ste P- However, no significant increase in mass-loss rate 
below ~24kK is identified from our results. 

Finally, recall empirical values of the Lamers et al. (1995) 
'bistability jump' were adopted by Vink et al. (1999, 2000) 



in their radiatively driven wind calculations. These predictions 
have subsequently been used in evolutionary calculations (e.g. 
Meynet & Maeder 2000) and spectral synthesis calculations 
(e.g. Rix et al. 2004). Consequently, our re-determination of 
physical parameters and wind properties of early B supergiants 
has potential consequences for evolutionary and spectral syn- 
thesis calculations. 



5.4. UV morphology 

The primary role of the present study is to provide physi- 
cal parameters and wind properties of Galactic B supergiants 
based on optical diagnostics, with reference to the bistability 
'jump'. We shall postpone a detailed discussion of UV spec- 
tral comparisons to elsewhere, but wish to comment on one 
specific aspect relating to UV morphology. Walborn (1971) in- 
troduced the B0. 7 subtype on the basis of the behaviour of Si m 
/14552/Siiv ,14089. Walborn & Nichols-Bohlin (1987) found 
that B0. 5 la and B0. 7 la supergiants were well defined and 
quite distinct in their UV spectral morphology, with regard to 
the presence/absence of Cn and Aim P Cygni lines, plus the 
shape of the Si iv absorptions. Between B0. 5 and B0. 7 the blue- 
shifted Si iv absorptions suddenly become narrower and deeper, 
the latter reversing the smooth trend from earlier types. 

We have shown that the bistability 'jump' does not have a 
well defined edge, but rather a general trend to lower Vco/vesc, 
particularly between B0. 5 and B0. 7 subtypes. How successful 
are our models, which are based on optical diagnostics, at re- 
producing these distinct UV properties? 

In Figs. |9jab) we present comparisons between de- 
reddened IUE (HIRES short wavelength SWP) observations of 
selected B supergiants and synthetic spectra. The UV spectra of 
B supergiants are greatly affected by metal blanketing. Indeed, 
-95% of the strong (W A > 0.01A) lines in the /LU250-1900 re- 
gion of early B supergiants are due to Fe m-rv transitions. Fe iv 
dominates the far-UV blanketing at B0, equal numbers of Fe m 
and Fe iv lines contribute at B0. 7, whilst Fe in dominates at later 
subtypes. The density of spectral lines is so large that there are 
no clean UV continuum regions, such that our comparisons are 
between fluxed datasets. 

In general, the UV spectra are fairly well reproduced by 
the optically derived diagnostics, both prominent P Cygni wind 
profiles of Si iv, C iv and Al m, and the primarily Fe-blanketed 
continuum. The most notable exceptions are (i) Civ /L11548- 
51 for subtypes later than B0. 5 for which negligible P Cygni 
emission is predicted, in dramatic contrast with observation - 
N v /U1238-42 suffers a similar problem for all subtypes - and 
(ii) The P Cygni absorption strength of Siiv /LU393-1402 is 
predicted too strong in B0-0.5 supergiants. 

Observationally, C iv and C n lines co-exist in B0.7-B2 su- 
pergiants, whilst the the former is effectively absent for stellar 
temperatures below ~25kK. One obvious means of maintain- 
ing high ionization states at relatively low stellar temperatures 
(and hence Civ, Nv wind lines) is via X-rays, produced by 
shocks in the stellar wind, which are also responsible for the 
clumped winds discussed above. Quantitatively addressing this 
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Fig. 9. (a) Comparison between de-reddened HIRES IUE spectroscopy (/LU250-1575) of selected early B supergiants (solid 
black) and synthetic UV spectra resulting from optical analyses (red in electronic version, dotted in paper version). 
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Fig. 9. b) Comparison between de-reddened HIRES IUE spectroscopy (/L11575-1900) of selected early B supergiants (solid 
black) and synthetic UV spectra resulting from optical analyses (red in electronic version, dotted in paper version). 
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question is beyond the scope of the present paper, and will be 
investigated in a future study. 

Nevertheless, we also have a partial answer to the reversal 
of the conventional UV morphologies for two early B super- 
giants, namely the BO. 5 la supergiant HD 1 15842 and BO. 7 la 
supergiant HD 91943 (Walborn et al. 1995). In general, B0.5 
and BO. 7 supergiants are distinct in terms of Cn and Aim 
due to the substantial differences in ionization between ~26kK 
(B0.5) and ~22.5kK (B0.7) of these elements. Indeed, the un- 
usually high stellar temperature of ~24.5kK for HD 91943 
(BO. 7 la) implies a high degree of ionization, such that negligi- 
ble wind features are predicted, or observed, at C n and/or Al m. 
This high temperature results from an unusually low wind den- 
sity, relative to other BO. 7 supergiants, as evidenced by its weak 
Ha signature in Fig. 01. In contrast, HD 115842 (B0.5Ia) has 
a rather strong wind, and so a low stellar temperature, with the 
Al m feature seen weakly. 

Finally, we turn to the other morphological difference be- 
tween BO. 5 and BO. 7 supergiants in the UV, i.e. narrower and 
deeper Si iv absorption at later subtypes. We have demonstrated 
above that there is indeed a trend for stars with T e ff > 24kK - 
approximately dividing BO. 5 from BO. 7 subtypes - to possess 
higher Voo/v esc ratios, i.e. broad, shallow Si iv absorption is ob- 
served in BO. 5 stars with high wind velocities, with narrow, 
deep Si iv absorption observed in BO. 7 stars with low wind ve- 
locities. 

Omitting HD 91943 and HD 115842 from the sample of 
Howarth et al. (1997), the average BO. 5 and BO. 7 la supergiant 
wind velocities are 1535 km/s and 955 km/s, respectively. The 
wind velocities of the morphologically unusual stars are indeed 
unusual. HD 91943 (B0.7Ia) has a fast wind of 1470 km/s, 
whilst HD 1 15842 (B0.5 la) has a slow wind of 1 180 km/s. Of 
course, a physical explanation is incomplete since the present 
models fail to predict the observed behaviour of Si iv in cases 
with fast winds amongst BO-BO. 5 supergiants, i.e. all except 
HD 1 15842 show too strong Si iv absorption. 

Clumping may play a role here too. If the conventional 
mass-loss diagnostics (e.g. Ha) in B supergiants formed in re- 
gions where the wind is clumped, the global mass-loss rate 
will decrease, potentially desaturating Si iv absorption. Similar 
conclusions were reached by Evans et al. (2004b) in their far- 
UV study of late O and early B supergiants in the Magellanic 
Clouds, in which the observed Siv /LU062-73 absorption 
could only be matched by clumped models (see also Prinja et 
al. 2004). A test calculation has been carried out for HD 38771 
in which clumping was included, with a volume filling factor 
of 0.1. Si iv /H393-1402 was not found to de-saturate in this 
case (C iv /H548-51 did de-saturate), indicating either a higher 
clumping factor or an alternative explanation. 

6. Conclusions 

We present detailed optical studies of a sample of Galactic 
B supergiants in order to determine physical parameters and 
wind properties based on contemporary model atmosphere 
techniques, for comparison with recent results on Magellanic 
Cloud B supergiants (e.g. Evans et al. 2004b, Trundle et al. 
2004). With respect to Kudritzki et al. (1999) who carried out a 



similar study based on an adopted temperature scale and a Ha 
mass-loss determination from an unblanketed code, we obtain 
lower stellar temperatures (and hence luminosities) and either 
comparable (early B) or higher (mid B) mass-loss rates. CNO 
elements are found to be partially CNO processed in general. 
On average, [N/C] and [N/O] are increased relative to (recently 
revised) Solar abundances by +1.1 dex and +0.8 dex, respec- 
tively. The morphologically unusual B supergiant HD 2905 
(BC0.7 la) indeed has an unusually low N enrichment and C/O 
depletion. 

With respect to recent theoretical predictions, we find rea- 
sonable agreement with those of Vink et al. (1999, 2000) for 
early (>23kK) B supergiants assuming homogeneous winds, 
although the predicted winds for later subtype B are too strong. 
This would only be exacerbated in the case of clumped winds in 
the Ha line forming region (Evans et al. 2004b). With respect to 
recent SMC B supergiant studies, Galactic mid-B supergiants 
are observed to have systematically stronger winds, although 
the situation is less clear for early B supergiants. We obtain 
reasonable agreement between the mid-IR continuum derived 
mass-loss rates of Barlow & Cohen (1977) and Ha mass-loss 
rates suggesting similar clumping factors for these regions, in 
contrast with expectation (Runacres & Owocki 2002). 

We have investigated the so-called 'bistability jump' 
amongst early B supergiants, in which Lamers et al. (1995) 
claimed a step in the ratio Voo/ve S c at 21kK (Bl) from 2.65 to 
1 .3. From our detailed study of larger sample close to this sub- 
type, we find that there is in fact a gradual decline in v^/v^c 
from early to mid B supergiants (also identified by Evans et al. 
2004a for SMC supergiants) from -3.4 above 24kK, -2.5 for 
20-24kK, and -1.9 below 20kK. 

We compare our optically derived spectral fits to UV spec- 
trophotometry of early B supergiants. Morphological differ- 
ences in the UV between B0.5 (~26kK) and B0.7 (~23kK) 
subtypes (Walborn et al. 1995) are attributed to the difference 
in ionization and wind velocity between these subtypes. Recall 
Vink et al. (1999) predicted a shift in ionization around 25kK 
for B supergiants, due to which wind densities increase dramat- 
ically below this temperature. The former does indeed match 
the observed shift in wind velocities at ~24kK, although no sig- 
nificant increase in mass-loss rate is inferred from our results, 
affecting evolutionary models which make use of such predic- 
tions (e.g. Meynet & Maeder 2003). Finally, we discuss po- 
tential solutions to spectroscopic discrepancies at Si iv (at B0- 
0.5) and Civ (at B0. 7-1. 5) which again suggest the presence 
of shocked, clumped winds in the line forming regions. Further 
progress requires the development of spectroscopic clumping 
diagnostics. 
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